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ABSTRACT 

We discuss the relationship between a standard IShakura fc SunvaevI (|1973[ ) accretion 
disk model and the Big Blue Bump (BBB) observed in Type 1 AGN, and propose a 
new method to estimate black hole masses. Wc apply this method to a sample of 23 
radio-loud Narrow-line Seyfert 1 (RL-NLSl) galaxies, using data from WISE (Wide- 
field Infrared Survey Explorer), SDSS (Sloan Digital Sky Survey) and GALEX. Our 
black hole mass estimates are at least a factor above previous results based on 
single epoch virial methods, while the Eddington ratios are correspondingly lower. 
Hence, the black hole masses of RL-NLSl galaxies are typically above IO^Mq, in 
agreement with the typical black hole mass of blazars. 
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1 INTRODUCTION 

00 ' The Spectral Energy Distribution (SED) of Active Galac- 
tic Nuclei (AGN) spans several orders of magnitude in fre- 
quency and results from the superposition of radiation emit- 

■ ted by different components. 
In radio-quiet Type 1 sources, characterized by the 

^N) , presence of broad emission lines in their optical spectrum, 
the most luminous components are the "Big Blue Bump" 
(BBB, between ~1 /^m and ~3 nm, or log(!//Hz) ~14.5- 
- , , . 17) and the "infrared bump" (IR bump, between ~1 mm 

■ and ^1 jun or log(i^/Hz) ~11.5-14.5). The former is the 
5— ( \ most prominent feature in the SED (| Sanders et al.l 1 19891 : 
^ ■ lElvis et all 1 1994 [Richards et"alll200d) . while the latter ac- 
counts for 20-40% of the bolometric AGN luminosity. The 
BBB is thought to be thermal radiation from the accre- 
tion disk, while the IR bump is thermal radiation emitted 
from a dusty torus located a ~1 pc from the black hole 
l|Sanders et al.|[l989l ). Superimposed to the BBB there is of- 
ten a minor component named "Small Blue Bump" (SBB, 
extending from 2200 A to 4000 A) which is likely the blend- 
ing of several iron lines and hydrogen Balm er continuum 
I Wills et all 1 19851 : IVanden Berk etall bOOll V This scheme 
roughly describes the SED of AGN o ver at least 5 orders 
of magnitude in bolometric luminosity (|Sanders et al.|[l989l : 
lElvis et all Il994l : iRichards et al.] bOOd ). It also applies to 
powerful blazars, although in these cases two more com- 
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ponents are needed to describe the entire SED: the "syn- 
chrotron hump" (extending from radio to IR/optical wave- 
lengths) and the "Compton hump" (extending from X- 
rays to TeV energies) which may overwhelm the torus and 
the BBB radiation. These further components characterize 
radio-loud sources whose jet is closely aligned to the line of 
sight, and are due to the synchrotron and inverse Compton 
processes, respectively. 



The common energy production process in AGN is be- 
lieved to be accretion onto a super-massive black hole (Af ~ 
10^~ ^"Mq), through a disk whose observational properties 
depend (among other parameters) on the black hole mass 
and accretion rate. This interpretation led several authors 
to use mo dels of geometrically thin, optically thick accre- 
tion disks (|Shakura fc Sunvaevlll973l . hereafter AD model) 
to fit the optical/UV SED of AGN in order to determine the 
black hole mass and the accretion rate (e.g 
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Malka n fc SargcntI Il982l : iMalkanI Il983l : IZheng et al 
Sun fc Malkan 198i). The AD fitting method allowed to 
estimate such quantities for those active nuclei which are 
too distant {z > 0.1) or too bright for other direct meth- 
ods, such as resolved s tellar/gas dynamics, to be applied 
(Pcrrarcsc & Ford"2005V_ However, as discussed extensively 
in Koratkar & Blacs (.19991 ) and references therein, such sim- 
ple models provide only rough fits to the observed data. 
Among the major issues with this interpretation, we point 
out a few ones: the broad-band continuum slopes a,^ (with 
Fi, oc v" ") at optical-NUV wavelengths found in litera- 
ture fe.g. iNeugebauer et al.lll979l : IVanden Berk et al.ll200ll : 
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iDavis et al] l2007l : iBonning etHI 120071 ') are incompatible 
with the slope = 1/3, expected from the AD model; the 
spectrum from a simple accretion disk does not reproduce 
the observed power law extending at X-rays and the soft 
X-ray excess (e.g. iPounds et al.lTl98^ : iNandra fc Pounds! 
1 1994 iFabian fc Miniuttill2005l ): the gross properties of the 
spectrum o f radio-quiet AGN appear to scale wit h the 
luminosity jSanders et al.lll989l: IWalter fc Fiiii3ll993l '). but 
does not shift in frequency ( Laor fc David 1201 la ). The lat- 
ter issue indicates that the BBB spectrum peaks always 
at, or near, the same frequency. However, this observa- 
tion is hard to reconcile with reasonably broad distribu- 
tions of black hole masses, accretion rates, inclinations and 
radiative efficiencies. Given these difficulties, the AD fit- 
ting method is not widely employed as a black hole mass 
estimation method. Rather, it is sometimes used to in- 
directly infer other parameters such as the accr etion ef- 
ficiencv jDavis fc LaoJ 1201 ll : iLaor fc Da^ l2011al . but see 
iRaimundo et al.ll2012l')" or to explain specific quasar proper- 
ties (|Laor fc Davisll20l"lbl '). while the black hole mass is usu- 
ally estimated using reverberation mapping calibrated scal- 
ing relatio ns (the so-called "single epoch viria l method" , or 
SE y i rial. |Petersonlll993l: IPeterson et al.ll2004l: lOnken et al.1 



12004 IVestergaard fc Petersonll2006l : iBentz et al.ll2009l l. 

The SE virial method is currently emplo yed to esti- 
mate black hole masses in large catalogs (e.g. IShen et al.l 
I2OIII ) for its simple applicability. However, the resulting 
black hole mass estimates are known to be affec ted by un- 
certainties of a factor of ~3 (jPark et al.l I20l3 ). In addi- 
tion, there may be systematic uncertainties related to the 
Broad Line Region (BLR) geometry and inclination; the 
role of radiation pressure; the modeling of emission line 
profiles; the contribution from other components (e.g. host 
galaxy); intrinsic differences between different AGN; and 
whether gravity dominates the motions of the BLR clouds 
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Marconi et "al]|2008l : iDecarh et al]|2008l . I2OIII : 



[Peterson 2011 ). Black hole masses estimated in this way are 
subjected to a number of assumptions and may therefore be 
just order of mag nitude estimates (e.g. ICroo m 2011. but see 
lAssef et al.ll2012l ). 

There has been claims about the existence of a corre- 
lation between t he SE virial mass and th e radio luminosity, 
either absolute (|Franceschini et alii 19981 ) or normalized to 
the optica l luminosity (i.e. the radio loudness parameter, 
lLaodl2000l ). It was further noticed that a black hole mass 
greater than ~ 3 x 10* Mq would be required in order to 
develop a relativistic jet as observed in powerful blazars. 
Later, these findings have been rev i sited w ith the availabil- 
ity of larger samples. IWoo fc Urrvl (|2002l ) found that both 
the radio-quiet and radio-loud AGN span the same range of 
black hole mass, and that there is no evidence for a strong 
correlation between the radio-loudness and black hole mass. 
Recently, th e issue on the black hole mass threshold has been 
revisited bv lChiaberge fc Marconil (|201ll ): by taking into ac- 
count the radiation pressure on the BLR clouds in comput- 
ing the SE virial masses they found that a black hole mass 
> 10* is required to produce a radio-loud AGN. 

In this paper we revisit the AD spectrum fitting method 
and show (^ that the AD model provides a rather satis- 
factory description of the Type 1 AGN SED in the major- 
ity of cases (at least at optical/NUV wavelengths) once the 
contributions from other emitting components (such as host 



galaxy and/or jet) have been properly taken into account. 
Therefore, the AD modeling method is a viable and indepen- 
dent way to estimate black hole masses. This is particularly 
interesting for a class of AGN sources for which black hole 
masses are suspected to be systematically underestimated 
by the S E virial method: the Narrow Line Seyfert 1 (NLSl) 
galax ies (|Marconi et al.ll2008l : iDecarh et alll2008l : IPetersonI 
I20T1I) . 

NLSl sources are characterized by the relatively small 
values of the full width at half maximum (FWHM) of the 
"broad" component of the H/J emission line (FWHM(H/3) < 
2000 km s"^), by the presence of stro ng blended iron lines, 
and of a promine nt soft X-ray excess (|Osterbrock fc Poggel 
ll985l : |Poggj|2000l ). By estimating the virial black hole mass 
using the H/3 emission line, and [O in] wi dth as a surrogate 
for th e bul ge stellar velocity dispe rsion, iGrupe fc Mathun 
(|2004 ) and (|Mathur fc Grupell2005l ) claimed that NLSl lie 
systematically below the M-a* of non-active and active 
Broad Line galaxies (BLSl). This indicates that the black 
hole masses of NLSl are systematically smaller than the 
black hole masses of BLSl for a given value of u*. The same 
considerations apply when considering objects of the same 
luminosity: NLSl appear to accrete at a higher Eddington 
ratio with respect to B LSl, with some ob jects exceeding the 
Eddington luminosity (|Zhou et al.ll2006l ) . 

Recently, a few NLSl sources have been confirmed to 
be part of a ne w class of 7-ray emi tt ing sources, as detected 
by Fermi/hAT (|Abdo et al.ll2009al l3: lFoschinill201ll ). besides 
blazars and radio-galax ies. Variability of the 7-ray emission 
(|Calderone et al.]|201ll ) allows to exclude a starburst origin 
of the 7-rays and confirms the presence of powerful relativis- 
tic jets, such as those found in typical blazars. The emerging 
picture is that 7-NLSI sources are very similar to powerful 
blazars, except for their small widths of broad lines, and con- 
sequently their sm all SE virial black hole masses (10®~*Mq, 
I Yuan et al.ll200^ ). Hence these sources are the best candi- 
dates to settle the question on whether very massive black 
holes (> 10* Mq) power relativistic jets. If the SE virial 
mass estimates will be confirmed this would imply that a 
large mass is not required to produce a radio-loud AGN. 
On the other hand, if masses turn out to be systematically 
under-estimated in NLSl, we would then conclude that ex- 
tragalactic jetted sources preferentially live in large black 
hole mass systems. 

Furthermore we may find that NLSl too lie on the 
"canonical" M-cr* relation of broad line sources and accrete 
close to (albeit below) the Eddington rate. Possible expla- 
nations for the observed small widths of the broad lines in 
NLSl include: the virial mass scaling relations may need to 
be modified in order to account for radiation press ure effects 
(IMarconi et al. I I2OO8I : IChiaberge fc Marconill201lf ): the BLR 
may have a disk-like geometry and be oriented almost face- 
on, so that the Doppler shifted line velo city, projected alon g 
the line of sight, turns out to be small (iDecarli et al.ll2008l ): 



a combination of these effects (|Peterson 2011 



In this work we show that the broad-band composite 
SEDs of AGN are roughly compatible with a simple, non- 
relativistic AD model (^Ijl, and discuss a method to estimate 
the total disk luminosity using either the continuum ( ^2.1|l 
or the line luminosities ( H2.2p as proxy. Then, we perform a 
spectroscopic analysis of the SDSS spectra of 23 radio-loud 
NLSl (UJ in order to disentangle the host galaxy and/or jet 
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contribution from the AGN continuum, and estimate line 
luminosities f i]3.ip . In ^ we discuss a new method to es- 
timate the black hole mass and accretion rate, using AD 
spectrum modeling, and apply this method on the afore 
mentioned sample. We discuss our results in ij5]and|6l and 
dra w our conclusion in j} 71 Th e observational properties of 
the IShakura fc SunvaevI (|l973l ) AD model are summarized 



in the appendix. 

Throughout the paper, we assume a ACDM cosmology 
with Ho = 71 km s"^ Mpc"\ = 0.27, = 0.73. 



1.1 Notation 

In what follows we will consider a non-relativisticQ steady 
state, geometrically thin, optically thick accretion disk 
IShakura fc Sunvaevlll973 '). extending from i?in = 6iig to 
Rout = 2 X lO'^-Rg, where 7?g — GM/c? is the gravitational 
radius of the black hole. The integrated disk luminosity is 
Ld = / L^dv = riMc?, with 77 ~ 0.1 (radiative efficiency). 
The corresponding Eddington ratio is ^ = id/i/Edd with 
I/Edd = 1.3 X 1O'''^(M/1O^M0) erg s'V 

The relation between the disk luminosity and its 
"isotropic equivalent" counterpart is Vl° = {2cos9)Ld (Eq. 
IA13[) . where 6 is the angle between the normal to the disk 
and the fine of sight. For Type 1 AGN we take (2 cos 6) = 1.7 
(Ea.[KW\i. 

We will refer to the peak frequency in the vL^ repre- 
sentation as i/p, and to the luminosity of the peak as VpL^-p- 
These quantities scale with the physical parameters as fol- 
lows (EQ.[Mland lX9| : 

Vp on M ' M ' 
i^pL^p oc M. 



In particular we notice that u^Lvp ~ 0.5Ld (Eg. lAlO)) . The 
location of the peak (i.e. its luminosity and frequency) de- 
termines uniquely the black hole mass and accretion rate. 
Details about the observational properties of the AD spec- 
trum are given in appendix lAll 

All spectral slopes oiv are defined as F,j oc v"" . 



2 ACCRETION DISK SPECTRUM IN AGN 
SPECTRA 



[Richards et~al] (|2006l ') built an average Type 1 QSO SED 
using data from 259 (mainly radio-quiet) sources, observed 
with instruments ranging from radio wavelengths to X-rays. 
Individual SEDs have been interpolated between available 
bands. An average SED is then computed as a geometric 
mean of individual ones, and is shown in Fig. [1] (red line). 
Also shown in F i g. [Jar e: a spiral galaxy template as given in 
iMannucci et al.l (|200ll . orange), normalized to have a bolo- 
metric luminosity of 10"^^'^ erg s~ ^ ; the l ocation of the Small 
Blue Bump (SBB, IWills et all Il985l : IVanden Berk et al.1 
I2OOII ): three reference frequencies corresponding to 5100A, 
3000 A and 1350 A (red filled circles), commonly used in cal- 
culation of bolometric luminosity (see below); the average 



^ General relativistic correction are negligible for the purpose of 
our work, see ^A4I 



spectral slopes found in literature, as measured on compos- 
ite spectra at near IR, optical/UV (*Vandcn Berk et~al ]|200ll . 
green) and far UV wavelengths ((Tclfcr ct al. 2002|. purple); 
the rest frame frequency range covered by SDSS, for values 
of z =0, 0.3, 1, 2 and 3 (thin blue lines). Finally, we show 
the AD spectrum that best fits the composite Type 1 SED 
at optical/UV wavelengths (black line). The parameter for 
the AD model are: log(M/M0)=9, ^=0.05 and 6) = 30°. 

The agreement between the AD model and the compos- 
ite SED is rather good, therefore the association between the 
BBB and thermal emission from simple AD model is justi- 
fied, at least in the interval 1000-5000A, or log(///Hz) = 14.8- 
15.5 (black dotted vertical lines). A few discrepancies be- 
tween the AD model and the composite Type 1 QSO SED 
arise: 

• at log(://Hz) <14.7 a further component emerges in the 
spectrum, which may be either the host galaxy, the emission 
from a dusty torus or some other component; 

• at log(i'/Hz) ~15 a Small Blue Bump (SBB) is present, 
likely due to a blending of iron lines and Hydrogen Balmer 
continuum; 

• at log(!//Hz) >15.6 other physical components con- 
tribute to the fiux (e.g. a corona). 

Note that, in this interpretation of the BBB, the portion of 
the AD spectrum characterized by the = 1/3 slope (thick 
blue line) is hidden by the host galaxy and the torus com- 
ponents, and cannot be revealed directly with observations 
(although in some case it may be detected in polarized light, 
iKishimoto et al.ll2008l ). The average slopes at optical/UV 
and far UV (green and purple lines) are roughly consistent 
with the slopes near the peak of the AD spectrum. We no- 
tice however that fixed spectral features (such as the SBB) 
may affect the estimation of spectral slopes. Furthermore, 
the value of the slope likely depends on the width of the 
wavelength range inside which it is defined. Therefore it is 
not always possible to infer the presence of an AD spectrum 
by just checking the spectral slopes at optical/UV wave- 
lengths. By contrast, at near IR the average slope (green 
line) is inconsistent with an AD spectrum, but this is likely 
due to the host galaxy component. 

We conclude that the AD model provides a reasonable 
description of the gross properties of Type 1 AGN SED at 
optical/NUV wavelengths, and the similarity between the 
predicted spectrum and the average BBB is rather strong. 
Under this assumption it is possible to infer the black hole 
mass and the accretion rate by comparing the observed SED 
with the AD spectrum, as discussed in SJ?) Our black hole 
mass estimation method requires an estimate of the disk 
luminosity {Ld), as discussed in the following two sections. 



2.1 Continuum luminosity as a proxy to disk 
luminosity 

The broad-band similarity among AGN spectra allows to 
use the continuum luminosity at a given wavelength as a 
proxy for the bolometric luminosity, that is Lboi = Cboi x 
\L\. In order to explore this relationship [Richards et al.l 
(|2006h measured the bolometric luminosity for each spec- 
trum (defined to be the integral isotropic luminosity between 
lOO^m and 10 keV) and derived a bolometric correction 
(Cboi) based on the continuum luminosity at 3000A, 5100A 
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AD model: M=10^ M©, Edd. ratlo = 0.05, = 30^ 
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Figure 1. Comparison of the composite Type 1 AGN SED (red line) from [Richards et al.l ||200(J ) an d an AD model fill. Ill for 
log(M/MQ)=9, £=0.05 and d = 30° (black line). Also shown are: a spiral galaxy template as given in iMannucci et al.l ll200 l|. or- 
ange), normalized to hav e a bolometric luminosity of 10*^'^ erg s"'^; the location of the Small Blue Bump (SBB, Iwiuset al.lll985l : 
I Vanden Berk et al.ll20o"lh : three reference frequencies corresponding to 5100A, 3000A and 1350A (red filled circles), commonly used 
in calculatio n of bolometric luminosi ty; the average spectral slopes f ound in liter ature, as measured on composite spectra at near IR, 
optical/UV l lVanden Berk et al.ll20o"ll . green) and far UV wavelengths l lTelfer et al.l 2002. purple); the rest frame frequency range covered 
by SDSS, for values of z =0, 0.3, 1, 2 and 3 (thin blue lines). Thick blue line highlights the portion of the AD spectrum characterized 
by the slope = 1/3. The rest frame frequency range inside which the AD model reproduces the shape of the AGN composite SED 
(log(!^/Hz)=14.8— 15.5) is shown with dotted black lines. 



and 3/im. The resulting values are: Cboi(3000A) = 5.62 ± 
1.14, Cboi(5100A) = 10.33 ± 2.08, Cboi(3/im) = 9.12 ± 2.62. 
The distribution of Cboi values is relatively narrow, with 
a relative dispersion of the order of ~20%. Note however, 
that in particular cases the Cbo i value can d iffer b y as much 
as 50% from the mean value. IShen et al.l (|201ll . hereafter 
Sll catalog) have shghtly re-calibrated the Cboi values, and 
extended the analysis to 1350A, in order to compute bolo- 
metric luminosities for all the sources in their sample. Their 
values are: Cboi(5100A) = 9.26, Cboi(3000A) = 5.15 and 
Cboi(1350A) = 3.81. 

In order to calibrate analogous relations to estimate the 
disk luminosity Ul° we numerically estimate the bolometri c 
luminosity of the composite SED in [Richards et all (|2006l ). 
and compare it with the disk luminosity for the AD model 
shown in Fig. [T] The resulting relation is: 

-^-d" ~ g-^'bol (1) 



Then, we compare Ul° with the luminosities at 5100A, 
3000A and 1350A wavelengths, as measured on the com- 
posite SED: 

LT ~ 4.4 i/L^(5100A) 

LT ~ 2.4 iyL,(3000A) (2) 

LT ~ 1.8 i/L,(1350A) 

The locations of these wavelengths are shown with red filled 
circles in Fig. [1] Considering the uncertainties (~20%) of 
Cboi we conclude that our relations (Eq. [T] and [2| are com- 
patible with those of Sll. Eq. [5] provides a reliable estimate 
of LT as long as the source continuum is not dominated 
by other emitting components such as host galaxy starlight 
or synchrotron radiation from a relativistic jet (for radio- 
loud sources). In these cases we need alternative luminosity 
estimators, as discussed in the following section. 
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2.2 Line luminosities as a proxy to disk luminosity 

Relations similar to Eq.[2]can be obtained by using line lumi- 
nosities. Line ratios ar e known to be approximately constant 
among AGN ((Francis ot al, 1991. : Vand en Berk et al. 20o3): 
by setting the Lya luminosity to 100, relative luminosities 
of H/3, Mgll and Civ (both narrow and broad components) 
lines are 22, 34 an d 63, respectively, while the total line lu - 
minosity is 555.8 (|Francis et al.lll99ll : ICelotti et all 1 19971 ). 
Therefore it is possible to have a rough estimate of the lu- 
minosity of all emission lines by measuring the luminosity of 
a single line. Also, according to the photo-ionization model, 
the line-emitting gas is ionized by the accretion disk con- 
tinuum radiation. Therefore we expect (to a first approxi- 
mation) the disk to line luminosity ratio to be a constant: 
Ul° = KLiino. This provides a way to estimate the disk lu- 
minosity using a single (or a few) line luminosity estimates. 
In order to calibrate the k parameter we consider all sources 
in the Sll catalog having both a continuum and line lumi- 
nosity estimate for at least one of the combinations: 5100A- 
H^, 3000A-Mgll and 1350A-Civ. The number of sources 
in each subsample are 22644, 85514 and 52157 respectively 
(note that a single source typically belongs to two such sub- 
samples). For each source we estimate Luno using the broad 
and nar row line luminosities give n in Sll and the coeffi cients 
given Inl^andseTaD (119911 ) and lCelotti et all ([l993). Then 
we compute the k parameter as follows: 

_ i^r(Eqia> 

K — 

where the disk luminosity Ul° is computed using the contin- 
uum luminosity given in Sll, and Eq.[2l The distributions of 
ft for the three combinations are approximately log-normal 
(Fig. [21 upper panel) with median values: 

log «;(5100A - H/3) = 1.08 ± 0.28 

log/t(3000A - Mgii) = 1.10 ± 0.21 (4) 

log k(1350A - C iv) = 0.92 ± 0.28 

The widths of the k distributions in Fig.[2]show that the disk 
luminosity LJf° computed using the continuum and the line 
intensities differs by <0.3 dex, i.e. a factor <2. Hence, the 
relationships between continuum and line luminosities seem 
quite robust. A possible explanation for the larger disper- 
sion in the 5100A-H/3 case, with respect to the 3000A-Mgn 
one, may be that the continuum luminosity at 5100A is con- 
taminated by the host galaxy. 

Both the continuum and the line luminosities in Sll are 
affected by uncertainties, therefore the distributions shown 
in Fig. [5] are likely broadened by measurement errors. Thus, 
the intrinsic dispersion is expected to be smaller than 0.3 
dex (a factor of ~2) for H/3 and Civ, and 0.2 dex (a factor 
of ~1.6) for Mgll. This is yet another evidence that SEDs in 
most AGN show some degree of universality: the constancy 
of the continuum to line luminosity ratio at optical wave- 
lengths implies a constant optical continuum to ionizing UV 
luminosity ratio. 

Since the samples in the Sll catalog are dominated by 
radio-quiet sources (the great majority are undetected in 
the FIRST survey), we repeat our analysis on radio-loud 
sources, i.e. those sources for which the radio-loudness pa- 
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Figure 2. Upper panel: distribution of the k parameter (Eq. |3] 
for the three combinations 5100A-H/3, 3000A-Mgn and 1350A- 
C IV. Both the continuum and line luminosity estimates are those 
reported in the Sll catalog. The number of sources in each sub- 
sample are 22644, 85514 and 52157 respectively (note that a single 
source typically belongs to two such subsample). Lower panel: the 
same as upper panel, for the subsamp le of Narrow-Line Seyfert 1 
sources common to both the Sll and I Zhou et al. I 120061) catalog. 

ramete]0is greater than 100. Interestingly, the k parameters 
for the radio-loud sub-sample of Sll differ by at most ~5% 
from the values quoted above. 

By using the values given in Eq. [4] and the coefficients 
to compute the line luminosity Z/uno discussed above, we are 
able to estimate the total disk luminosity as follows: 

Lr = 12L(H/3)^ = 3031(H/3) 

Lr = 12.5L(Mgii)^= 204i(Mgii) (5) 

LT = 8.4L(Civ) = 74.1L(Civ). 

63 

We repeat the above analysis on the subsample of NLSl 
sources, that are the focus of our study . In particu l ar, we 
consider the sources common to both the lZhou et al.l (|2006l ) 
and the Sll catalog (1210 sources). The distribution^ of 
the K parameter are still log-normal, and are shown in the 
lower panel of Fig. [5] Median values are now ~0.15 dex (i.e. 
a factor ~1.4) greater: 

log k(5100A - H^) = 1.23 ± 0.23 
log k(3000A - Mg ll) = 1.24 ±0.21. 

The uncertainties are of the same order of magnitude. The 
class of NLSl sources is therefore characterized by both a 



^ The radio-loudness parameter provides an indication of 
whether the AGN SED is dominated by radiation at radio fre- 
quencies or optical band. Historically, it is defined as the rati o 
of 5 GHz to optical B-band luminosity l lKellermann et al]|l989l) . 
The values we used here are those given in Sll, defined as the 
ratio of flux densities at 6 cm and 2500A (rest frame). 
3 The 1350A-C IV case is missing since the SDSS wavelength cov- 
erage does not allow to observe both the C IV line and the H/3 line 
(required to classify the source as a NLSl). 
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smaller width and a smaller luminosity of lines. The result- 
ing disk luminosities are: 



We dropped spec- 
one of the following 



LT = 424L(H/3) 
1)1° = 286L(Mgii). 



(7) 



In order to estimate the accretion disk luminosity using 
a single spectrum we can use either Eq. [2] (whose uncertain- 
ties are ~20%) or Eq. [7] (whose uncertainties are a factor 
~2). In cases where the observed continuum radiation is 
dominated by components other than AD, e.g. synchrotron 
emission from the jet or host galaxy starlight, Eq. [2] would 
overestimate the disk luminosity. Therefore Eq. [7]is our pre- 
ferred choice to estimate Ul° ■ 



3 THE SAMPLE 

The aim of this work is to estimate the black hole masses of 
the sample of 23 Rad io-Loud, Na r row-L ine Seyfert 1 sources 
(RL-NLSl) given in I Yuan et all (|2008l . hereafter Y08). We 
identify each source with a sequential index (#1, #2, etc.), 
following the same order as in Tab. 1 of Y08. 

All sources have been spectroscopically observed in the 
SDSS, and 21 over 23 sources are also in the Sll catalog. 
The IR photometry at 3.4^m, 4.6/im, 11.6/jm, 22.1/jm from 
WISE (Wide-field Infrared Survey Explorer, IWright et al.l 
I2OIOI ) is available for all sources. Finally, 21 o ver 23 sources 
have photometric measurements by GALEX (|Martin et al.l 
I2OO5I ). either in the Medium Imaging Survey (MIS), or the 
All sky Imaging Survey (AIS). We noticed that when multi- 
ple GALEX observations were available, we find significant 
variability in a few cases (#3, #8, #18) possibly due to the 
jet component. In these cases we chose preferably the MIS 
photometry with lower luminosity. 

The redshifts are in the range z =0.1-0.8, therefore the 
continuum in the SDSS spectra will likely trace the AD com- 
ponent The FWHM(H;g) are less than 22 00 km s~\ 
as req uired by the definition of NLSl given in IZhou et al.l 
(|2006l ). The SE virial black hole masses are in the range 
log(M/M^) =6-8 , while the Eddington ratio are (. =0.5-3 
I Yuan et al.ll200§ ). The radio morphology is compact, unre- 
solve d on 5" scale, and the radio loudness (|Kellermann et ahl 
1 19891 ) is >100 for aU sources. 

The overall observational pro perties are very similar to 
that of blazars (|Yuan et al.ll200^ ). and the 7-ray emission 
from these sources has been predicted, and later detected 
in 7 RL-NLSl sour ces (|Abdo et al.ll2009d : ICalderone et"al] 
l201ll : lFosch"inill201ll '). 4 of which are in the Y08 sample. How- 
ever, these sources show unusually small widths of broad 
emission lines, and consequently small SE virial black hole 
mass estimates, when compared to typical blazars. 

In order to apply our black hole mass estimation method 
(311 to the sources in the sample we need to perform a spec- 
troscopic analysis of the SDSS data. In particular we need 
to disentangle the host galaxy and/or jet contribution from 
the AGN continuum, and estimate the emission line lumi- 
nosities. This procedure is described in the following section. 



3.1 Spectral analysis 

We used the s pectra from the Sloan Digital Sky 
Survey (SDSS, lYork et all I2OO0I ). data release 7 



(DR7, lAbazaiian et al.l [200i). 

tral bins marked by at least 

mask flagslfl SP_MASK_FULLRE JECT , SP_MASK_LOWFLAT , 
[SP_MASK_SCATLIGHT , SP_MASK_BRIGHTSKY , SP_MASK_M ODATAl 
'SP_MASK_CQMBINEREJ, .SP_M ASK_BADSKYCHI Also, we 
dropped 100 bins at the beginning and end of each spec- 
trum, in order to eventually avoid artifacts from instrument 
or pipeline. 

Each spectrum has been de-reddened using the Galac- 
tic extm£tion__ra estimated from dust IR emission maps 
in ISchlegel et al. (1998 ), an d the extinction law reported 



in lCardeUi et all (| 19891 ') and lO'DonneUl (|l994h . We are cur- 
rently neglecting any intrinsic reddening in the rest-frame 
of the source. Then we transformed the spectra to the rest 
frame by assuming isotropic emission (i.e. multiplying the 
flux by 47r_D£). The redshift estimates are provided by the 
SDSS pipeline. Finally, we rebinned each spectrum by a fac- 
tor of 3 in order to improve the signal to noise ratio, result- 
ing in a spectral resolution of \/5\ ~1450 (corresponding 
to ~200 km s"^). 

The model used to fit the spectra consists of five com- 
ponents: 

• a smoothly broken power law to account for the AGN 
continuum ("AGN continuum" component) ; 

• a spiraffl host galaxy template spectrum from 



[m annucci et al.1 (|200ll ') and a power law to (eventually) ac- 
count for the synchrotron emission from the jet ("galaxy" 
and "jet" components respectively). The galaxy compo- 
nent has a single free parameter (the overall normalization) 
The parameters for the jet component are estimated us- 
ing data from WIS E (Wide-field Infrared Survey Explorer_ 
I Wright et al.|[201ol ). In particular we use the photometry in 
the two bands at the longest wavelengths (ll^m and 22/im 
to estimate the luminosity and the slope of the power law 
If the resulting slope is greater than -1 we extrapolate the 
power law to optical wavelengths and subtract the contribu- 
tion from the SDSS spectrum. Otherwise we do not consider 
any jet component. Parameter of the jet component are fixed 
during the fitting process; 

• the iron templates from IVestergaard fc WilkesI (l200l[) 
(at UV wavelengths) and from Veron-Gettv et al. ( 20041 ') (at 
optical wavelengths); 

• a Gaussian profile for each emission line listed in Tab. 
[1] The FWHM of narrow lines are forced to be in the range 
200-1000 km s~^, while that of broad lines are forced in 
the range 1000-3000 km s~\ Furthermore, the FWHM and 
velocity offset of the H/3 narrow component is tied to the 
width and offset of [O ill] A4959 and [O ill] A5007. 

• a maximum of 10 additional Gaussian line profiles 
which are not "a priori" associated to any specific transi- 
tion. These components are necessary to account for (e.g.) 
the iron blended emission lines in the range 3IOO-35OOA (not 
covered by the above-cited iron templates) , or line asymme- 



* See |http: //www. sdss . org/dr7/dm/f latFlles/spSpec .html[ 

^ The results of the spectral fitting procedure do not change 
significantly by consid ering the elliptical galaxy template from 
iMannucci et al.l l l200lD . 

6 In analyzing the source SDSS J094857.32+00222 5.5 (#5) we 
also applied an exponential cutoff at logfi^/Hz] = 14 llAbdo et al.l 



I2OO9 



: 
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Table 1. List of emission lines used in modeling SDSS spectra. 
Third column (Type) indicates if a broad (B), a narrow (N) or 
both components are used in the fit. 



Line 



Wave [A] Type I Line Wave [A] Type 



C II 


2326 


B 


N I 


5199 


N 


Mg II 


2798 


BN 


He I 


5876 


BN 


Ne V 


3426 


N 


Fe VII 


6087 


N 


O II 


3727 


N 


O I 


6300 


N 


Ne III 


3869 


N 


Fe X 


6375 


N 


H 5 


4101 


B 


N II 


6548 


N 


H 7 


4340 


BN 


H a 


6563 


BN 


O III 


4363 


N 


N II 


6583 


N 


He II 


4686 


BN 


S II 


6716 


N 


H/3 


4861 


BN 


S II 


6731 


N 


O III 


4959 


N 


Ar III 


7136 


N 


O III 


5007 


N 









tries. The FWHM of the additional hnes are forced to be in 
the range 1000-3000 km s~^, except for hnes in the range 
3100-3500A for which the upper hmit is lO'* km s^^. A pos- 
teriori, we check whether the wavelength range identified by 
the full width at half maximum of these additional hnes con- 
tains any of the transition hnes listed in Tab.[T] In this case 
we associate the two components, and numerically compute 
the line luminosity on the composite line profile. 

Results of the spectral fitting are shown in Tab. [2] and 
Fig. Ill 



4 BLACK HOLE MASS ESTIMATION 
METHOD 

The AGN continuum in the rest frame wavelength range 
1000-5000A (or log(:/ /Hz) = 14.8-15.5), if interpre ted as ra- 
diation emitted from a lShakura fc SunvaevI (|l973D accretion 
disk (^, allows to constrain an AD model, and to infer the 
black hole mass. Once we assume proper values for the in- 
ner radius of the disk 7?in and the viewing angle 6 ( i]4.ip . 
the luminosity and frequency of the peak of the AD spec- 
trum uniquely identify a value of the black hole mass. In 
the following sections we will discuss two methods to locate 
the peak of the AD spectrum, and infer the black hole mass 
and accretion rate. An example of the application of both 
methods to a specific case will be discussed in t)4.3l 



4.1 Hypotheses 

The methods relies on the following hypotheses, which need 
to be independently verified: 

(i) accretion in AGN occurs through steady-state, geo- 
metrically thin, optically thick, non-relativistic accretion 
disks. The emitted spectrum is well described by an AD 
model (|Shakura fc SunvaevlflQT^ '): 

(ii) once the galaxy and/or jet contribution has 
been subtracted, the continuum radiation in the range 
log(i//IIz) = 14.8-15.5 (311) is emitted directly from the ac- 
cretion disk, i.e. it has not been reprocessed by intervening 
material, nor it is emitted by some other component; 

(iii) the spatial extention of the disk is -Rin = 6-Rg , corre- 
sponding to a radiative efficiency r] ~ 0.1. The outer radius 



of the disk -Rout = 2 x 10 -Rg is not critical, since at frequen- 
cies much smaller than the AD spectrum will always be 
hidden by other emitting components. The assumption for 
-Rin, on the other hand, is more critical, since our black hole 
mass estimates show a linear dependence on this value (case 
(v) of ilX2l) : 

(iv) the relation between disk luminosity and its 
"isotropic equivalent" counterpart is L'J° = (2cosS)L^ (Eq. 
IA16|I . Since we are interested in Type 1 AGN the viewing an- 
gle is i n the range 0-45 deg (i. e. the aperture of the obscuring 
torus. ICalderone et al.l |2012| '|. The averaged de-projection 
factor is thus (2cosS) ~ 1.7 (Eg. lAlSl IA16|) . corresponding 
to a viewing angle of ^^30 deg. 

The AD model has four parameters: M, M, -Rin and cos 9 
( i]A3|l . With the assumptions discussed above, the remaining 
unknown parameters are the black hole mass M and the 
accretion rate M. 



4.2 Procedure 

Usually the localization of the peak of the AD spectrum is 
not accessible by using a single instrument, requiring opti- 
cal/UV multiwavelength observations. When these observa- 
tions are available it is possible to constrain the AD model, 
and estimate the frequency Vp and luminosity VpL^^° of the 
peak. The latter can then be used to infer the total disk 
luminosity LJf" (Eg. lAlOl) . Finally, the black hole mass and 
the accretion rate can be estimated as follows: 



M 



W^Mq 

M 

Mq yr- 



1.44 



= 0.21 



1015 Hz 



r IS. 



(2 COS 6') X 10*5 erg s" 



1/2 



ris. 



(2 COS 6*) X 10« erg s- 



(8) 

The uncertainties on these results can be estimated by prop- 
agating the uncertainties in the and i^pL^^° parameters in 
the above equations. Hence, whenever the data allow to con- 
strain the location of the peak of the AD spectrum, the ac- 
curacy of the black hole mass estimate is determined only b y 
the accuracy of the data points (e.g. ISbarrato et al.l [20121 ). 
When UV observation are not available (or not reliable) the 
location of the peak cannot be constrained, and we must 
resort to an alternative method. 



The LINE 



Here we propose a new method for the AD modeling which 
relies on broad line luminosities to estimate the total disk 
luminosity. Fig.Oillustrates the method. We use Eq.[7lto es- 
timate -LJf°. When both the H/3 and Mgll line luminosities 
were provided by our spectral fitting ( t)3.ip we considered 
the average of the resulting disk luminosities. This enables 
us to estimate a value for the luminosity of the peak fpL^^" 
(Eg. lAlOll . i.e. to fix a "ceiling" in the z^L^ representation 
(black dashed line in Fig. [Sj: the peak of the AD spectrum 
must lie on this line. Then we use observations from a sin- 
gle instrument (SDSS) to constrain the peak frequency Vp, 
which is related to the black hole mass. In particular, we shift 
the AD spectrum horizontally (green arrow), until the AD 
spectrum reproduces the AGN continuum identified in ij3.1l 
Note that the model to be compared with data in Fig. Ois an 
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Table 2. Results of the spectral fitting for the 23 RL— NLSl sources in lYuan et al ] l|2008l ') catalog. Columns are: (1) source numeric 
identifier; (2) SDSS name of the source; (3) redshift; (4) luminosity and error of the H/3 emission line (both the broad and narrow 
components); (5) luminosity and error of the Mgll emission line (both the broad and narrow components); (6) wavelength Aq and (7) 
luminosity Xq^Xq used to constrain the LINE model (see i|4.2.1ll ; (8) jet component (extrapolated from WISE data to wavelength Aq) to 
AGN continuum luminosity ratio. 
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"isotropic equivalent" AD spectrum (Eg. IA13| . The result- 
ing AD model (red solid line) provides an estimate for !/p, 
to be used (along with -£/d°) in Eq. [8]to infer the black hole 
mass and accretion rate. Finally, we cornpute the Eddington 
ratio i using the luminosity of the disk|3| I = I/d/^Edd- 

The main source of uncertainty in the process is the un- 
certainty in the disk luminosity (a factor 2, ^2.2[) . This 
uncertainty is shown as a grey shade in Fig. [S] In order to 
evaluate a confidence interval for our estimates of M we re- 
peat the whole process requiring the peak luminosity of the 
AD spectrum to lie respectively at the top and the bottom 
of the grey stripe. The resulting AD models (dot-dashed 
red lines) provide respectively the lower and upper limits of 
the confidence interval on the black hole mass, which typi- 
cally is ±0.5 dex. In some case these limiting AD models are 
too distant from the data to provide a meaningful descrip- 
tion of the AGN continuum. This occurs typically for the 
low luminosity solution, corresponding to the upper limit in 
black hole mass (e.g. i]4.3p . In these cases a visual inspection 
would reduce the thickness of the grey stripe, and hence the 
uncertainty on the black hole mass. 

Further sources of uncertainties are the assumption on 
the radiative efficiency 77 ~ 0.1 and on the viewing angle 
S ~ 30 deg ( H4.ip . The uncertainty due to the former can be 



^ We are neglecting the contribution from the torus in computing 
the Eddington luminosity since it is reprocessed radiation from 
the disk. If we had used the bolometric luminosity Lboli instead 
of -L(ji the resulting value would be overestimated by a factor 3.4 
on average (Eq.lTland lATet . 



estimated by considering that our black hole mass estimate 
is M oc 77 (case (v) of i]A2|) , and that the actual value of 77 is 
expected to range from ~6% (for non-rotating black hole) to 
at most ~30% (for a spin parameter a = Jc/GM^ = 0.998, 
lThornelll974l ). Therefore the uncertainty on the black hole 
mass due to the uncertainty on 77 (and ultimately on the 
black hole spin) is +0.5/-0.2 dex. If we consider the possi- 
bility that the black hole can be maximally counter-rotating 
(with respect to the direction of accretion) then the uncer- 
tainty on the black hole mass due to the uncertainty on 
77 becomes ±0.5 dex. The uncertainty due to the assump- 
tion on the viewing angle can be estimated by propagating 
the error in Eq. [S] Typically, this is negligible compared to 
the uncertainties discussed above, being at most 0.04 dex (a 
factor ~1.1), provided ^max <45 deg. The uncertainties due 
to Ul° and 77 are likely uncorrelated, therefore, the maxi- 
mum expected uncertainty for the black hole mass estimate 
is ±0.7 dex. 

The LINE procedure can be implemented without any 
fitting procedure, provided we have an estimate for the 
broad line luminosities and the AGN continuum ( tj3.ip . The 
search for the peak frequency can be implemented by iden- 
tifying a wavelength Ao and the corresponding luminosity 
of the AGN continuum Ao-£/Ao, and requiring the AD spec- 
trum to match this luminosity at the same wavelength. A 
comparison between the resulting AD model and the AGN 
continuum can be performed "a posteriori" in order to as- 
sess the reliability of the black hole mass estimate ( i]4.2.2p . 
For SDSS spectra of sources with z < 0.8 the value of Ao 
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Figure 3. The LINE black hole mass estimation procedure: the 
SDSS source spectrum (black line) is analyzed with our fitting 
procedure f ilS.ll l in order to estimate the broad line luminosities. 
Then we use Eq.[7]to estimate 1/^°. This is equivalent to estimate 
a value for the luminosity of the peak i/pL^° (Eg. lAlOt . i.e. to fix 
a "ceiling" in the vLi, representation (black dashed line). Then 
we use the SDSS spectrum to constrain the peak frequency Vp, 
which is related to the black hole mass. In particular, we shift the 
AD spectrum horizontally (green arrow), until the AD spectrum 
reproduces the AGN continuum identified in i|3.1l The resulting 
AD model (red solid line) provides an estimate for i/p, to be used 
(along with ijjf°) in Eq.|8]to infer the black hole mass and accre- 
tion rate. The uncertainty in the disk luminosity LJj*" (a factor 2, 
i|2.2ll is shown as a grey shade. In order to evaluate a confidence 
interval for our estimate of M we repeat the whole process requir- 
ing the peak luminosity of the AD spectrum to lie respectively 
at the top and the bottom of the grey stripe. The resulting AD 
models (dot-dashed red lines) provide respectively the lower and 
upper limits of the confidence interval on the black hole mass. 

has been chosen empirically as follows: 

Ao ~ Amin I -r 1 

\ Amin / 

where Amin and Amax are the minimum and maximum rest- 
frame wavelengths of the SDSS spectrum. This value is suf- 
ficiently close to the short wavelength edge in order to min- 
imize the contamination from other continuum components 
(either galaxy or jet); Ao is also sufficiently far from the 
shortest available wavelength, at which the estimated lumi- 
nosity may be unreliable due to noise and/or edge artifacts. 
Following these prescriptions the LINE method can be ef- 
ficiently implemented as an automated procedure on large 
samples. 

^.2.2 The BEST procedure 

In order to assess the reliability of the LINE procedure we 
proceed with a visual localization of the AD spectrum peak 
using the SDSS and GALEX observations. Photometry from 
GALEX has been de-reddened following the same proce- 
dure as for the optical SDSS spectra ( jjS.lf) . In addition, 
when a jet component is considered in the spectral fitting, 
we compute the jet-subtracted GALEX photometry. We can 
not exclude that further absorption took place either in the 
AGN environment or the intervening medium, therefore we 



consider the photometry as lower limits to the actual rest- 
frame luminosity. Note that SDSS and GALEX data are not 
simultaneous, therefore it may happen that these data sets 
trace the source in two dif ferent state of emission , e.g. a disk 
or a jet dominated state (jCalderone et al.]|2012h . 

For each source we require the slope of the AD spectrum 
to match as close as possible the slope in the AGN contin- 
uum (^STj, and to lie above the (jet-subtracted) GALEX 
photometry. A few exceptions to these rules will be consid- 
ered in ij6] Since this is a manually tuned AD model we give 
no error associated to the corresponding black hole mass. 

4.3 Example of application of the methods 

As an example we discuss the case of SDSS 
J09531.7.09+283601.5 (#6), in Fig. H The WISE pho- 
tometry is shown with black filled circles. The spectral fit 
( ^3.1|l is shown as a black line, while the AGN continuum 
component is shown as cyan line. The jet power law 
extrapolation from IR data is shown as a purple line. The 
GALEX photometry and their jet subtracted counterparts 
are shown as open circles and symbols respectively. 

The disk luminosity LJf° with its uncertainty of a factor 2 is 
shown as a grey stripe: the peak of the AD model must lie 
within this region. The grey dashed grid shows the location 
of peaks for AD models with values of black hole mass and 
Eddington ratio shown respectively below and above the 
grey stripe. 

By applying our black hole mass estimation methods 
we identify the LINE and BEST AD models, shown with a 
red and orange solid line respectively. Both AD models pro- 
vide a rather good representation of the AGN continuum. 
The BEST AD model, however, needs a slightly higher lumi- 
nosity than the LINE model in order to lie above the (jet- 
subtracted) GALEX photometry. Note that the observed 
spectrum (black line) has a significantly lower spectral slope 
(i.e. it is "redder") than the AGN continuum, because of the 
host galaxy and jet contributions. Having considered these 
components in the spectral analysis allows us to reveal the 
real AGN continuum (cyan line) whose slope agrees with our 
AD spectrum. 

In order to evaluate the uncertainty on the LINE black 
hole mass we repeat the procedure by requiring the AD 
model to peak at the top and the bottom of the grey stripe. 
The resulting AD models (shown with dot-dashed red lines) 
are found to bracket the real case: the lower one cannot ac- 
count for the AGN continuum, while the higher one is signif- 
icantly above the GALEX photometry. This situation often 
occurred during the analysis of the sources (^SJ, therefore 
our black hole mass uncertainties are rather conservative. 

Our AD models can be compared with those corre- 
sponding to the SE virial masses and bolometric luminosities 
reported in the Y08 and Sll catalogs (green and blue lines). 
We consider the disk luminosity as computed using Eq. [T] 
Note that our peak luminosities are very similar to those 
of Y08 and Sll, since this is the condition we required (on 
average) to calibrate Eq. [T] However, these models do not 
provide a good description of the AGN continuum because 
their peak frequencies lie ~0.25 dex above our estimates of 
^p, therefore our black hole mass estimates are 0.5 dex (a 
factor --^3, Ea. lA8|l greater than the virial ones. The possible 
reasons to explain such differences will be discussed in i]6.1l 
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Figure 4. Application of the black hole mass estimation methods to the source SDSS J09531. 7. 09+283601. 5 (#6). The WISE photometry 
is shown with black filled circles. The spectral fit ( ^3.11 is shown as a black line, while the AGN continuum component is shown as cyan line. 
The jet power law extrapolation from IR data is shown as a purple line. The GALEX photometry and their jet subtracted counterparts 
are shown as open circles and "+" symbols respectively. The disk luminosity with its uncertainty of a factor 2 is shown as a grey 
stripe: the peak of the AD model must lie within this region. The grey dashed grid shows the location of peaks for AD models with 
values of black hole mass and Eddington ratio shown respectively below and above the grey stripe. The LINE ( i|4.2.1ll and BEST ( i|4.2.2ll 
AD models are shown with a red and orange solid line respectively. In order to evaluate the uncertainty on the LINE black hole mass we 
repeat the procedure by requiring the AD model to peak at the top and the bottom of the grey stripe. The resulting AD models are 
shown with dot-dashed red lines. 



5 RESULTS 

We analyzed the data from the 23 RL-NLSl sources of the 
Y08 catalog. The spectral analysis (^STTJ of each individual 
source is shown in Fig. IBll The results are summarized in 
Tab. [5] The fitting models are in good agreement with data 
with reduced in the range 1.16-1.86. Also, the jet contri- 
bution at optical/NUV wavelengths is typically negligible, 
except for the #2, #11, #15, #16, #20 and #22 sources. 

The results of our black hole mass estimation methods 
are shown graphically in Fig. ICll (adopting the same 
notation as in Fig. 3)). The results are summarized in Tab. 
[3] The AD models identified by the LINE procedure provide 
a rather good description of the AGN continuum in 17 over 
23 cases (indicated with a blank in the second column of 
Tab. [3|. The remaining 6 sources cannot be modeled with 
an AD spectrum, and are considered "bad cases" (indicated 



with a "*" symbol in the second column of Tab. |3] see ^ 
for a discussion of these sources). These sources will not be 
considered in the following analysis. 

The comparison between the black hole mass estimates 
for the LINE and BEST AD models are shown in Fig. (5] The 
mean value for the ratio of the two mass estimates is: 



log 



M[LINE] 
A/ [BEST] 



0.07 ±0.37 



(9) 



The two black hole mass estimates are therefore compatible, 
within the uncertainties associated to the LINE procedure 

(mj. 

In Fig. [6] we show the comparison between the black 
hole masses from the AD models (LINE in upper panels, 
BEST in lower panels) and the black hole masses from SE 
virial method, as given in the Y08 (left panels) and Sll 
(right panels) catalogs. The uncertainty associated to SE 
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Table 3. Results of our black hole mass estimation method. Columns are: (1) source numeric identifier; (2) flag to indicate if the AD 
"signature" (i.e. the slope > —1 a t optical wavelengths, see Q is missing. (3) peak frequency of the AD model, (4) black hole mass 
estimate (with its uncertainties) and (5) Eddington ratio for the AD model identified by our automatic procedure (LINE model); (6), (7), 
(8) corresponding quantities for the BEST model; single epoch (SE) virial black hole mass estimate given in the (9) Y08 and (10) Sll 
catalogs. 
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virial mass is assumed to be 0.5 dex. In Fig. [7] we show the 
histogram of the ratio of our black hole mass estimates to 
the SE virial ones from the Y08 (left panel) and Sll (right 
panel) catalogs. The mean values for the ratio of the mass 
estimates are: 

/, MfLINEl \ /, AffLINEl \ 

/, MfBESTl \ /, A/ [best! \ 

(10) 

The mean values are of the same order (or even greater) 
than the maximum uncertainty associated to the LINE black 
hole mass estimate (0.7 dex, i ]4.2.1|) . therefore our black hole 
mass estimates are not compatible with the SE virial ones. 



6 DISCUSSION 

As discussed in ^ the characteristic disk spectral slope 
= 1/3 cannot be directly observed in AGN SED. However 
for values of the black hole mass log(Af/AfQ) > 8 and Ed- 
dington ratio £ < 1, the peak of the AD spectrum is (in prin- 
ciple) observable. Indeed, for 17 over 23 sources considered 
here, the SDSS continuum show an increasing trend in the 
fL,^ representation (slope > —1) at log(i//Hz) > 14.8, 
where the accretion disk is expected to dominate over other 
emitting components. Assuming that the observed BBB is 
actually radiation emitted directly from the accretion disk, 




7.5 8.0 8.5 9.0 

log(M/MQ) [BEST] 



9.5 



Figure 5. Comparison between black hole mass estimates ob- 
tained using the LINE and BEST procedures. Since the BEST esti- 
mate is a manually tuned AD model we give no error associated 
to the corresponding black hole mass. 
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the slope > —1 at log(i^/Hz) > 14.8 becomes the "sig- 
nature" of the presence of the AD component. In the con- 
sidered sample 17 sources over 23 show such signature 
(the remaining six "bad" cases will be discussed below). Our 
spectral fitting procedure f i|3.1[l reveals an emission compo- 
nent that is well described by an AD model, although this 
has not been included "a priori" in our fitting model. There- 
fore, observational data are in agreement with hypothesis 
(i), as discussed in t|4.1l Furthermore, the absorption col- 
umn densities A'h, as estimated from X-r ay spectral fitting , 
are compatible wi th the Galactic values (|Yuan et al.ll200a : 
iGrupe et al]|201(]| '). Hence, we expect the radiation we ob- 
serve not to be re-processed by any intervening medium, as 
assumed in hypothesis (ii) (*J4T]). Furthermore, we expect 
the contribution from other emitting components, such as 
host galaxy or jet, to be negligible at frequencies where our 
AD models are constrained (log(!//Hz) > 14.8). In particu- 
lar, the slopes in the AGN continuum component (qi/ > — 1) 
is inco mpatible with the ones i nferred from the galaxy tem- 
plate of iMannucci et al] |20o3). Also, the jet component is 
expected to decay at frequencies above a cutoff frequency of 
log(i//IIz) < 15, as in typical Flat Spectr um radio-quasars 
or powerful blazars (|Ghisellini et al.ll201Cll V 

The AD models identified by the LINE procedure 
f ^4.2.1l) for the 17 "good" sources provide a rather good 
description of the AGN continuum (Fig. IClj ). In particular, 
the AGN continuum slopes in the frequency range covered 
by SDSS are in good agreement with the ones from LINE 
AD models (red solid lines). Also, the two limiting solutions 
(dot-dashed red lines) likely bracket the real case, providing 
a robust estimate of our uncertainties. The average uncer- 
tainty on the black hole mass estimates are of the order of 
±0.5 dex (Tab. [3]). By taking into account the uncertainties 
due to hypotheses (iii) and (iv) (radiative efficiency 77 ~ 0.1 
and viewing angle ^ ~ 30 deg) we obtain a maximum un- 
certainty of ~0.7 dex ( {14.2.111 . 

In order to further assess the reliability of the LINE black 
hole mass estimates we considered the BEST AD models, 
identified by visually tuning the Ul° and Up parameters in 
order to achieve the best possible match between the AGN 
continuum identified in i|3.1l and the GALEX photometry. In 
a few cases we had to relax these requirements, as discussed 
below: 

• the assumption that we can reliably estimate the jet 
contribution at optical/UV wavelengths by extrapolating a 
power law from the WISE photometry ( i]3.ip may not be 
correct. For the #15 and 7^16 sources this assumption does 
not apply since the power law extrapolation (purple line) lies 
above the WISE photometry at shorter wavelengths (note 
that the error bars are smaller than the symbol in the plot). 
Source #5 would also falls in this class if the cutoff of syn- 
chrotron radiation at ~ 10^* Hz ( i]3.ip is neglected, since jet 
extrapolation would intercept optical SDSS data. In order 
to identify the BEST AD model for the #15 and #16 sources 
we used the continuum observed in SDSS data (black solid 
lines in Fig. IC1[) rather than the jet-subtracted one (cyan 
lines) as requirement at optical wavelengths. For #16 we 
obtained a good agreement between the BEST model and 
GALEX photometry. Lack of such agreement for #15 will 
be discussed below. A similar situation (jet component over- 
estimated at optical/NUV wavelengths) possibly occurs also 



for source #22. In order to be conservative, for this source we 
retained the original constraints to identify the BEST model. 
For the other sources the jet extrapolation is marginal at 
optical/NUV wavelengths (Tab. [2|, hence the assumption 
discussed here has a negligible effect. 

• for the #5 and #15 sources the GALEX photome- 
try does not follow the extrapolation from the SDSS slope. 
Therefore a single AD model is not compatible with both the 
SDSS and GALEX observations. This may be a consequence 
of source variability, since the SDSS and GALEX data are 
not simultaneous. Indeed, we found significant variability 
in GALEX photometry in a few sources ( tj3.ip . In particu- 
lar, so ur ce #5 is known to b e a variable source (|Abdo et al.l 
l2009al lbl: [Foschini et aLlbOllh . In these cases the BEST model 
is computed relaxing the requirement of taking GALEX 
photometry into account, and using only the SDSS data as 
guidelines. 

• for the #8 and #21 sources the SDSS and GALEX data 
appear to trace the peak of the AD spectrum. For these 
sources we neglected the jet component. Note that the BEST 
AD model for these sources provide a robust estimate of the 
black hole mass, since the peak of the AD spectrum has been 
directly observed. 

The BEST AD models are in good agreement with the LINE 
ones. In particular, note that the peak in BEST AD spectra 
lie inside the grey stripes for all sources except #17. Hence, 
the Eq. [7] are well calibrated. The resulting BEST black hole 
mass estimates are compatible the LINE ones (Eq. Fig. 
[5]). Also, the scatter in Fig. [5] (0.4 dex) is compatible with 
the uncertainty on the LINE estimates due to our ignorance 
on Ul° (0.5 dex, t|4.2.1|) . This provides further support for 
the reliability of the LINE black hole mass estimates. We 
conclude that, under the assumptions discussed in ^4. H our 
LINE procedure provides a reliable estimate of the black hole 
mass, within the quoted uncertainties. 

In six cases the LINE method do not provide an accept- 
able description of data (sources marked with a "*" symbol 
in the second column of Tab. O For these sources the ob- 
served SDSS continuum does not show an increasing trend 
(in the vL^ representation ) at log(!^/Hz) > 14.8. In two 
cases the SDSS continuum are dominated by the jet and/or 
host galaxy emission (#2 and #11), and the AD spectrum is 
not directly visible. In four cases (#7, #9, #10 and #20) the 
observed SDSS continuum appears "flat" in the vLi, repre- 
sentation, with no hints for a change of slope. Although the 
jet-subtracted continuum (cyan line) suggests the presence 
of an AD spectrum, this decomposition strongly depends on 
the assumption that the extrapolation of the jet component 
from IR data is also valid at optical wavelengths. In order 
to be conservative, we mark these sources as "bad", and 
neglect them in our analysis. 

6.1 Comparison with SE virial mass estimates 

The comparison with the SE virial mass reveals a systematic 
discrepancy between our mass estimate and those from the 
Y08 and Sll catalogs (Fig. |6]and[7ll. Although the discrep- 
ancy (>0.7 dex, Eg. llOf) is of the same order of the maximum 
uncertainty associated to the LINE procedure ( H4.2.ip . it ap- 
pears systematic. Therefore our black hole mass estimates 
are not compatible with the virial ones given in the Y08 
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Figure 6. Comparison between black hole mass estimates of AD models identified by our procedures (LINE models, upper panels, BEST 
models lower panels) and single epoch (SE) virial masses as given in the Y08 (left panels) and Sll (right panels) catalogs. The uncertainty 
associated to SE virial mass is 0.5 dex. 



and Sll catalogs. On average, our black hole masses turn 
out to be a factor ~6 (0.8 dex) greater than virial ones. 
A possible explanation for this black hole mass discrepancy 
may involve a radiative efficiency which is a factor ~6 lower 
than assumed. However, a value of ~ 0.02 is lower than 
the minimum efficiency expected for accretion onto a maxi- 
mally counter-rotating black hole (~ 0.03). Notice that the 
AD spectrum suggests that the accretion disk is still in the 
"radiatively efficient" regime; half the gravitational energy 
gained by matter at each radius is locally emitted as radia- 
tion, i.e. it is not advected into the hole. 

The mass discrepancy is not due to an inaccurate esti- 
mation of the jet contribution at optical/NUV wavelengths. 
If the actual jet contribution is lower than estimated, the 
corresponding AGN continuum luminosity (AoLaq (Tab. [2]) 
would correspondingly be higher. The "ceiling" luminosity 
level, on the other hand, are not affected by the presence of 
the jet, since it relies on line luminosities. In order to repro- 
duce an higher AqLao i retaining the same peak luminosity, 
the LINE AD model must shift to lower frequencies. There- 
fore, if the actual jet contribution is lower than estimated, we 
would have obtained greater LINE black hole mass estimates, 
and greater discrepancy with SE virial masses. Furthermore, 
the mass discrepancy is not due to having neglected the gen- 



eral relativistic corrections in the AD model. As discussed in 
i]A4l the AD model used throughout this work mimics the 
more sophisticated general relativistic one with ?7gr ~ 0.1, 
as long as frequencies below the peak are concerned. 

We speculate that a possible explanation for the mass 
discrepancy is a selection effect in calibrating the SE virial 
method. The virial method relies on the calibration of both a 
BL R size-continuum luminosity relation and of a virial fac- 
tor IPark et al.ll2012h . However, the sample used to perform 
the calibration consists of a few dozens of sources: the ones 
that has been reverberation mapped. As a consequence, the 
calibration of the method may be biased by selection effects. 
In particular, the method may provide significantly under- 
estimated black hole masses if the BL R has a flat disk-lik e 
geometry, and it is seen almost face-on (jDecarli et al.|[201ll ) . 
If these conditions apply, then the discrepancy between our 
mass estimates and the virial ones would be greater (on aver- 
age) for AGN showing the smallest widths of broad em ission 
lines, i.e. the class of NLSl sources (|Decarli et al ]|2008l ). The 
black hole mass estimates provided by our method, on the 
other hand, are scarcely affected by the viewing angle and 
the geometry of the BLR ( H4.2.ip . 

The Eddington ratios are in the range I =0.04-0.2 (Tab. 
[3]), significantly below the values reported in Y08. This dis- 
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Figure 7. Histogram of our black hole mass estimates to the SE virial ones from the Y08 (left panel) and Sll (right panel) catalogs. 
Our estimates are greater than SE virial ones by ^^1 dex for Y08 and ~0.8 for Sll (Eg. llOt . 



crepancy is due both to our greater black hole mass esti- 
mates (a factor ~6) and to the fact that we used the disk 
luminosity Ld (instead of Lboi) to compute the Eddington 
ratio (a factor --^3.4, Eg. [T] and IA16|) . Hence, our values of 
Eddington ratio a factor ^ 20 smaller than SE virial ones 
on average. With such small values of I the r ole of radiation 
pressure in determining the SE vir ial masses (|Marconi et al.l 
I2OO8I : IChiaberge fc Marconil 120111 ) is expected to be small, 
not sufficient to explain the black hole mass discrepancy. 



6.1.1 The "temperature" argument 

If the assumptions discussed in i]4.1l applv. then all the inde- 
pendently estimated black hole masses (even the SE virial 
ones) should produce an AD spectrum compatible with the 
observed data. This provides a simple way to compare our 
results with those reported in Y08 and Sll. The black hole 
mass discrepancy arises because the peak frequencies of the 
LINE and BEST AD spectra are significantly lower then the 
peak frequencies of the Y08 and Sll ones, while the peak 
luminosities are compatible. If the Y08 and Sll AD models 
(blue and green lines in Fig. lClfl were the correct ones, there 
must be a physical process able to shift photon to lower fre- 
quencies (i.e. to lower "temperatures") in order to account 
for the observed data. However, such a process cannot exist 
on a thermodynamic basis, since black body spectra (which 
build up the AD model) already have the lowest tempera- 
ture corresponding to a given luminosity and emitting sur- 
face. Since the luminosities are the same, the only way to 
reduce the temperature is to increase the emitting surface, 
that is by increasing the black hole mass. 



7 CONCLUSIONS 

In this work we analyzed the relationship between the Big 
Blue Bump (BBB) observed in the SED of Type 1 AGN 
and a lShakura fc SunvaevI (|l973l ) accretion disk model (AD 
model). The characteristic disk spectral slope = 1/3 can- 
not be directly observed in the AGN SED because of the 
contributions from other emitting components such as the 
host galaxy, the torus or (for radio-loud sources) the jet 
(32)1. Once these contributions are taken into account, the 
observations are compatible with the presence of an emitting 
component which is well described by an AD model. In par- 
ticular, the peak of such component can be observed directly 
in the frequency range log(!//Hz) = 1 4.8-15.5. By comparing 
the average Type 1 AGN SED from [Richards et"al] (|2006l ) 
with the AD model we calibrate the relations to estimate the 
total disk luminosity using the continuum line luminosities 
(at I350A, 3OOOA and 5100A) as proxy ( >i2.1|l . Furthermore, 
by us ing the emission line templates from iFrancis et al.l 
(|l99lD , we calibrate analogous relations based on the line 
luminosities of H/3, Mg 11 and C iv ( tj2.2|l . The latter provide 
more reliable disk luminosity estimates when the continuum 
is not dominated by the AD spectrum. 

The interpretation of the BBB as being due to the ther- 
mal emission from an AD allows to link the AGN observed 
properties to the properties of the super massive black hole. 
In particular, the luminosity and frequency at the peak of 
the AD spectrum uniquely identifies a value for the black 
hole mass and the accretion rate O. However, the direct 
observation of the peak of the AD spectrum requires broad- 
band multiwavelength observations. In order to estimate the 
black hole mass also for sources with a poor observational 
coverage, we propose a new method which relies only on 
spectral observations at optical wavelengths. In particular, 
we constrain the luminosity of the peak by using the line 
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luminosities. Then we constrain the frequency of the peak 
by requiring the AD model continuum and slope to repro- 
duce the observed AGN continuum beneath the emission 
lines ( i)4.2.1|) . The maximum uncertainty on our black hole 
mass estimates is ~ 0.7 dex (on average). This uncertainty 
is greatly reduced if the disk luminosity can be accurately 
determined, namely when the peak of the AD spectrum is 
visible within the frequency range of the data. 

We applied our method to the sample of 23 radio-loud 
Narro w line Seyfert 1 galaxies (RL-NLSl) of lYuan et al.l 
l|2008h . The method provides reliable black hole mass es- 
timates for 17 sources over 23, if our interpretation of the 
BBB is correct (fJS] O. The remaining six sources are ei- 
ther dominated by synchrotron radiation from the jet, or do 
not show "hints" for the presence of an AD-like emitting 
component. The resulting black hole mass estimates are a 
factor ~ 6 (on average) greater than the corresponding sin- 
gle epoch virial mass estimates, and the Eddington ratios 
are a factor ~ 20 below. The discrepancy between our black 
hole mass estimates and the single epoch virial ones may 
be due to selection effects occurred in the calibration of the 
BLR-continuum relation and of the virial factor. 

The black hole masses estimated in this paper for the 
sample of RL-NLSl are in the interval log(M/MQ) = 8-9, 
and the Eddington ratios are I =0.04-0.2. Therefore, very 
radio-loud NLSl appear not to be extreme in terms of black 
hole masses and Eddington ratios. Their black hole masses 
are similar to those of blazars. We find no evidence for jetted 
sources with mass below 10* Mq. 
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APPENDIX A: ACCRETION DISK MODEL 

Al Shakura&Sunyaev Accretion Disc (AD) 

Here we review the properties of the AD model 
for a geometrically thin, optically thick accretion disk 
jShakura fc Sunvaevlll973 . hereafter SS73 model) adopted 
in our analysis of the SED. 

The amount of gravitational energy released from each 
annulus of the disk is given by 
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where R is the distance from the black hole, Rg — GM/(? 
is the gravitational radius of the black hole, i?in is the inner 
radius of the disk. By introducing the adimensional param- 
eters: 
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we rewrite the emitted flux as F{R) = F{x) V, with 
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where all physical quantities are cast into the V parameter. 



while F{x) accounts for dimensionless flux distribution. The 
total disk luminosity is given by 

/-Rout 
2-kR F{R) dR 



-l/2x 



dx 



(A4) 



TjMc 



where Xout = RovLt/Rin is the normalized outer radius of 
the disk. The quantity in squared parentheses is equal to 1 
(provided -Rout 2> -Rin), therefore the parameter t] as deflned 
above is the radiative efficiency of the disk. By assuming 
-Rin = 6-Rg (appropriate for a non-rotating black hole) we 
obtain 77 ~ 0.1. 

The maximum amount of energy flux is (by differenti- 
ating Ea. lA3|) : 



MAX[-F(-R)] = F(-Rn 



64 

TT 



V 



(A5) 



and it is emitted at a radius -Rmax = 49/36 -Rin. The as- 
sumption of optical thickness implies that each annulus 
emits radiation as a black body with temperature: T(-R) = 
[F(-R)/cr]^/'*. The maximum temperature is therefore (Eq. 
[3: 



[K] 



3.46x10** 



0.1 



8/4 



M 



IO^Mq 



-1/2 



M 



Mq yr- 



1/4 



(A6) 

The emitted spectrum is a superposition of black body 
spectra: 

"-Rout 

U ^2x I 2tt RdR 7vB[u, T{R)] 

1/41 (A7) 



4Tv^RfnV 



,3/4 



X dx B 



Fix) 



where B[i>,T{R)] is the Planck function. The spectrum pro- 
flle is completely determined by the dimensionless integral, 
the only dependences on physical parameters [V) being the 
characteristic frequency (oc "p^/*) and the overall normaliza- 
tion (oc -Rfjj'P'^''''). The disk spectra are therefore self-similar, 
and the peak frequency and luminosity scale as: 



[Hz 



M 



-1/2 



M 



[erg s 



109 A/q 



M, 



yr" 



1/4 



M 



Mq yr- 



(A8) 



where is the frequency of the peak in the vLi, rep- 
resentation, log^ — 15.25 and logB = 45.36. By intro- 
ducing the Eddington ratio £ — Ld/LEdd (with LEdd = 
1.3 X 10'*^(-M/10^AjrQ) erg s~*), the previous equations can 
be rewritten as: 



[Hz] Vo.l 



L/2 / M 



V 109 Mq 



-1/4 



0.04 



1/4 



[erg s~^ 



M 
IO'JA/q 



0.04 



(A9) 



Notice that, for a given value of 77, an estimate of the lumi- 
nosity and peak frequency allows to determine the physical 
parameters M and M. 

The spectrum of an AD is shown in Fig. [T] (black solid 
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line): the superposition of black body spectra, weighted by 
the surface of emitting annuli produces a "flat" spectrum 
with slopes Qi/ ~ 1/3; at highest frequencies, the Wien spec- 
trum from the inner ring dominates, and the overall spec- 
trum decays exponentially. 

The self-similarity of AD spectra implies the existence 
of relations among quantities at the peak frequencies in the 
Lv and the vLu representations [v-pi and respectively), 
and the disk luminosity L^j: 



= 3.1 



= 0.66 



id 



= 0.5 (AlO) 



Also, note that the peak luminosity v-gLi,^ is independent 
from the actual value of -Rout, as long as Tiout > lO-Rin. The 
relation between the maximum temperature in the disk and 
the color temperature of the AD spectrum (i.e. the black 
body temperature associated to the peak frequency v^i) is 



3.5 Tcoi. 



A2 Peak shift 

The physical parameters 77, M and M uniquely identify the 
frequency and luminosity of the spectral peak (Eq. IA8I or 
IA9[I . Variations of one or more of these parameters will shift 
the peak along specific directions, whose slope in a log vL^ 
vs. log V plot is given by: 



^ _ dlogi/pL^p 
d log Vp 

Here is a list of peak shift relations used in this work 



(All) 



(i) vertical shift (q = oo): variations in M, M with con- 
stant MjM^ ratio (fixed r?); 

(ii) horizontal shift (a = 0): variations in Kl with con- 
stant M (fixed r\)\ 

(iii) a = 4: variations in A/, with constant M (fixed 77); 

(iv) a = —4: variations in M and M, with constant i oc 
M jM (fixed r;); 

(v) no shift: variations in all parameters, with constant 
TjM and ry/M. 

In particular, case (v) is used in H6.ll to show that if the 
actual radiative efficiency 77 is greater than hypothesized in 
i|4.1l then our black hole mass estimate is a lower limit. The 
physical interpretation of case (v) is depicted in Fig. lAll 
By following the black hole mass estimation method out- 
lined H4.2.1I we identify an AD model (black line) in agree- 
ment with observed data (brown line): the peak lies at the 
"ceiling" luminosity level determined by broad line luminosi- 
ties (Eq. [T]) and the AD spectrum is in agreement with the 
observed continuum. This model relies on the assumption 
of radiative efficiency 77 ~ 0.1. However, the actual value 
of efficiency may be different. By increasing the r\ param- 
eter (i.e. decreasing the inner radius of the disk Jilin) the 
peak shifts to higher frequencies and luminosities as radia- 
tion comes from the inner, hotter radii (blue line). This new 
model would no longer be in agreement with the "ceiling" 



° The slope oli, = 1/3 is achieved only for i?out — >■ oa. For a finite 
value of the outer radius of the disk, such as the one used in this 




log V 



Figure Al. The black hole mass estimate provided by the AD 
modeling procedures (Q show a linear dependence on the as- 
sumed value of the radiative efficiency: M ocrj (case (v) in ^A2| |. 
The physical interpretation of this dependence is shown in the fig- 
ure: the AD model (black line) identified with either the LINE or 
BEST procedure is in agreement with observed data (brown line). 
In particular the peak lies at the "ceiling" luminosity level deter- 
mined by broad line luminosities (Eq. [7]l and the AD spectrum is 
in agreement with the AGN continuum. By increasing the rj pa- 
rameter (i.e. decreasing the inner radius of the disk -Rin) the peak 
shifts to higher frequencies and luminosities as radiation comes 
from the inner, hotter radii (blue line). This new model would no 
longer be in agreement with the "ceiling" luminosity argument, 
therefore we must decrease the accretion rate (green line), leav- 
ing M and rj unchanged. Still, the obtained spectrum is not in 
agreement with the observed continuum, therefore we must de- 
crease the "temperature" of the spectrum ( N6.1.1I I. by increasing 
Al (red line). The final AD model is again in agreement with ob- 
served data, but has higher values of 77 and M (and a lower value 
of M). 



luminosity argument, therefore we must decrease the accre- 
tion rate (green line), leaving M and 77 unchanged. Still, the 
obtained spectrum is not in agreement with the observed 
continuum, therefore we must decrease the "temperature" 
of the spectrum f ^O.l.lfl . by increasing M (red line). The 
final AD model is again in agreement with observed data, 
but has higher values of 77 and Af (and a lower value of M). 



A3 Observational properties 

The emission from the whole (geometrically thin) disk is 
anisotropic since the observed flux is proportional to the 
projected area seen by the observer, i.e. oc cos 6', where 
is the viewing angle. By requiring f^^^ F,yD'ldQ, = we 
obtain: 



(1 -I- z) cos 61 



(A12) 



work -Rout 
av ~ 1/4. 



2 X lO'^-Rg, a more realistic value for the slope is 



where -Dl is the luminosity distance, Vo = v/il+z) is the ob- 
served frequency and F^^ is the observed flux density. Note 
that the luminosity-flux relation for a thin disk is different 
from the isotropic case, in particular the relation between 
the "isotropic equivalent" luminosity and the real luminos- 
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ity is: 

Lr = 2cos6lL^ (A13) 
The observed flux is therefore (from Eq. IA7I IA12|l : 

R dR 



rv„ ~ „ ( 1 + Z) cos 6 1 ■ ,, . 



(A14) 



The model for the observed spectrum has four parameters: 
M, M, Rin and cos 6 (the value of iiout is not important 
here) which are related to quantities in Ea. lA14l throua:h the 
temperature distribution given in Eq. lAll Not all parameter 
can be constrained observationally, since the viewing angle 
is degenerate with both Af and M. Hence we are forced to 
make a simplifying assumption about the inclination angle: 
since we are interested in Type 1 AGN, we assume that the 
viewing angle is in the range 0-45 deg, i.e. the aperture of the 
obscuring torus. If observed at a greater angle, the source 
would likely be classified as a Type 2 AGN. The average 
value of cos 9 (where 9 is measured from the disk normal) 
is: 



(cos 9) 



1 + cos 9n 



(A15) 



By setting Smax=45 deg we obtain (cos S)=0.854, corre- 
sponding to an average viewing angle of ~30 deg. With this 
assumption Eo. IA131 reads: 



1.7 U 



(A16) 



A4 General relativistic corrections 



The gener al relativistic model for th e accr etion disk is de- 
scribed in iNovikov fc Thornd (| 19731 ') and IPage fc Thornj 
1 19741 ). The differences with respect to the AD model in- 
fluencing the observational appearance of the spectrum are: 

(i) the innermost stable circular orbit (isco) depends on 
the spin parameter a — Jc/GM^. For a non-rotating black 
hole (a = 0) -Risco ~ 6R s ■ The maxim um spin of an accreting 
black hole is a = 0.998 (|Thorne|[l973 ). with Risco = 1.24i?g. 
The binding energy of a particle at fliaco i n units of the 
particle rest-mass is (e.g. ICunninghamlll975l l: 



2 Rq- 



3Ri, 



(A17) 



i.e. 77gr(a = 0) ~ 0.06 and r;gr(a = 0.998) ~ 0.32. This is 
expected to be the maximum possible value for the radiative 
efficiency (compare Eq. IA3[I . 

(ii) the different radial distribution of energy flux 
l|Page fc Thornd 1 19741 : IZhang et all \W9^ with respect to 
Eq. lAll The resulting spectrum is still the superposition of 
black body spectra; 

(iii) the spectrum received by distant observers is influ- 
enced by gravitational redshift, Doppler b oost and gravita- 
tional bending of light (|Cunningham|[l975h . 

iLi et al.l (|2005l ) have developed a package to synthesize the 
observed spectrum for an optically thick, geometrically thin 
accretion disk around a Kerr black hole, taking into account 
all these effects. The code is available as the mode l KERRBB 
within the X-ray data reduction package XSPEC (|Arnaudl 
Il996l ). In the following we will compare the spectral proflles 
of both the "classical" and "relativistic" models, and show 



that the differences are negligible for the purpose of our 
work. 

We compute the accretion disk flux, as received by an 
observer at a given distance, using both the the SS73 and 
KERRBB models. The black hole mass, accretion rate and dis- 
tance of the observer will be kept fixed for all the considered 
models. 

We consider five SS73 AD models, by varying the inner 
radius of the disk Rin- The values of Rin has been chosen in 
order to reduce the discrepancies between SS73 and KERRBB 
models (see below). We take the model with _Rin — 6-Rg as 
a reference spectrum, and normalize all other SS73 spectra 
by the luminosity of its peak {v^Fv ^^^). These spectra are 
shown with solid lines in Fig. IA2I Note that the only de- 
pendence on the viewing angle 9 for the SS73 model is due 
to the projected area seen by the observer, i.e. to a factor 
cos 6'. By plotting spectra normalized by cos 9 we completely 
remove this dependence. 

Then we consider five groups of KERRBB models, by vary- 
ing the spin of the black hole: a = —1,0,0.4,0.7 and 0.998. 
These values span the entire rang e of allowed values for 
the spin of an accreting black hole (|Thornelll974l ). For each 
value of the spin, we consider three different viewing angles: 
9 = Q° , 30° and 45° . All spectra are normalized by the lumi- 
nosity fp-Fy ^j,f of the reference SS73 model discussed above, 
and by cos9. These spectra are shown with dotted, dashed 
and dot-dashed lines in Fig. IA2I Note that the KERRBB mod- 
els show a residual dependence on the viewing angle, due to 
light bending and Doppler boosting. 

The values of 7?in has been chosen in order to allow 
the SS73 spectra to resemble as close as possible the KERRBB 
spectra, at given values of spin. The profile of the normalized 
spectra are indeed very similar (spectra of the same color), 
the differences being at most ±0.1 dex for the highest value 
of spin (a = 0.998). The (empirical) relation between the 
7?in in the SS73 model and the radiative efficiency of the 
corresponding KERRBB model is: 



Riri 

Rs 



2% 



+ 1.25 



(A18) 



where ijgr is given by Eq. IA17I 

From the observational point of view, the SS73 AD 
model with _Rin = 67?g (used throughout this work) mim- 
ics the KERRBB model with spin a ~ 0.7 {rjgi 0.1), as long 
as frequencies below the peak are concerned. Therefore, the 
results of our work are not influenced by having neglected 
the general relativistic corrections in modeling the accretion 
disk spectrum. 



REFERENCES 

Abazajian K. N., et al., 2009, ApJS, 182, 543 

Abdo A. A., et al., 2009a, ApJ, 699, 976 

Abdo A. A., et al., 2009b, ApJ, 707, 727 

Abdo A. A., et al., 2009c, ApJ, 707, L142 

Arnaud K. A., 1996, in Jacoby G. H., Barnes J., eds. Astro- 
nomical Data Analysis Software and Systems V Vol. 101 
of Astronomical Society of the Paciflc Conference Series, 
XSPEC: The First Ten Years, p. 17 

Asset R. J., Frank S., Grier C. J., Kochanek C. S., Denney 
K. D., Peterson B. M., 2012, ApJ, 753, L2 



© 2011 RAS, MNRAS OOO.[Tlf20l 



18 G. Calderone et al. 




-1.5 -1.0 -0.5 0.0 0.5 1.0 

log ^ / i^p(ref.) 

Figure A2. Comparison between the SS73 ijShakura fc Sunvaevll 197l3 . solid lines) and the KERRBB l|Li et al .120051. dotted, dashed and dot- 
dashed lines) accretion disk spectra. The SS73 spectrum with R^^ = 6 -Rg is used as a reference spectrum: all other spectra are normalized 
by its peak luminosity. We consider three viewing angles 9 = 0°, 30° and 45°, and normalize all spectra by cos0. This completely removes 
the dependence of the SS73 on the viewing angle. The KERRBB show a residual dependence on the viewing angle (dotted, dashed and 
dot-dashed lines). The colors identify a value for the inner radius of the SS73 model (respectively Ri^/Rg = 14.5, 10, 7.9, 6 and 2.8) and 
for the black hole spin of the KERRBB model (respectively a = —1, 0, 0.4, 0.7, 0.998). The inner radii for the SS73 models have been chosen 
in order to allow the SS73 spectra to resemble as close as possible the KERRBB spectra, at given values of the spin. The resulting empirical 
relation between iJin and the radiative efficiency of the KERRBB model (Eq. IA17|l is given in Eg. lAlsl 



Bentz M. C, Peterson B. M., Netzer H., Pogge R. W., 
Vestergaard M., 2009, ApJ, 697, 160 

Bonning E. W., Cheng L., Shields G. A., Salviander S., 
Gebhardt K., 2007, ApJ, 659, 211 

Calderone C, Foschini L., Ghisellini G., Colpi M., 
Maraschi L., Tavecchio F., Decarii R., Tagliaferri G., 2011, 
MNRAS, 413, 2365 

Calderone G., Ghisellini G., Colpi M., Dotti M., 2012, MN- 
RAS, 424, 3081 

Calderone G., Sbarrato T., Ghisellini G., 2012, MNRAS, 
p. L482 

Cardelh J. A., Clayton G. C, Mathis J. S., 1989, ApJ, 345, 
245 

Celotti A., Padovani P., Ghisellini G., 1997, MNRAS, 286, 
415 

Chiaberge M., Marconi A., 2011, MNRAS, 416, 917 
Groom S. M., 2011, ApJ, 736, 161 
Cunningham C. T., 1975, ApJ, 202, 788 
Davis S. W., Laor A., 2011, ApJ, 728, 98 



Davis S. W., Woo J.-H., Blaes O. M., 2007, ApJ, 668, 682 
Decarh R., Dotti M., Fontana M., Haardt F., 2008, MN- 
RAS, 386, L15 
Decarh R., Dotti M., Treves A., 2011, MNRAS, 413, 39 
Elvis M., et al., 1994, ApJS, 95, 1 

Fabian A. C, Miniutti G., 2005, ArXiv Astrophysics e- 
prints 

Ferrarese L., Ford H., 2005, Space Sci. Rev., 116, 523 
Foschini L., 2011 in Proceedings of "Narrow-Line 
Seyfert 1 Galaxies and their Place in the Universe". 
(arXiv: 1105. 0772) 
Foschini L., et al., 2011, MNRAS, 413, 1671 
Franceschini A., Vercellone S., Fabian A. C, 1998, MN- 
RAS, 297, 817 

Francis P. J., Hewett P. C, Fohz C. B., Chaffee F. H., 
Weymann R. J., Morris S. L., 1991, ApJ, 373, 465 

Ghisellini G., Tavecchio F., Foschini L., Ghirlanda G., 
Maraschi L., Celotti A., 2010, MNRAS, 402, 497 

Grupe D., Komossa S., Leighly K. M., Page K. L., 2010, 



© 2011 RAS, MNRAS OOO.[Tlf20l 



Mass of RL-NLSl 19 



ApJS, 187, 64 
Grupe D., Mathur S., 2004, ApJ, 606, L41 
Kellermann K., Sramek R., Schmidt M., ShafTer D. B., 

Green R., 1989, AJ, 98, 1195 
Kishimoto M., Antonucci R., Blaes O., Lawrence A., Bois- 

son C., Albrecht M., Leipski C., 2008, Nature, 454, 492 
Koratkar A., Blaes O., 1999, PASP, 111, 1 
Krolik J. H., 2001, ApJ, 551, 72 
Laor A., 2000, ApJ, 543, Llll 
Laor A., Davis S., 2011a, arXiv : 1110 . 0653 
Laor A., Davis S. W., 2011b, MNRAS, 417, 681 
Li L.-X., Zimmerman E. R., Narayan R., McCIintock J. E., 

2005, ApJS, 157, 335 
Malkan M. A., 1983, ApJ, 268, 582 
Malkan M. A., Sargent W. L. W., 1982, ApJ, 254, 22 
Mannucci F., Basile P., Poggianti B. M., Cimatti A., Daddi 

E., Pozzetti L., Vanzi L., 2001, MNRAS, 326, 745 
Marconi A., Axon D. J., Maiolino R., Nagao T., Pastorini 

G., Pietrini P., Robinson A., Torricelli G., 2008, ApJ, 678, 

693 

Martin D. C, et al., 2005, ApJ, 619, LI 
Mathur S., Grupe D., 2005, ApJ, 633, 688 
Nandra K., Pounds K. A., 1994, MNRAS, 268, 405 
Neugebauer G., Oke J. B., Beckhn E. E., Matthews K., 

1979, ApJ, 230, 79 
Novikov I. D., Thorne K. S., 1973, in Dewitt C., Dewitt 

B. S., eds, Black Holes (Les Astres Occlus) Astrophysics 

of black holes., pp 343-450 
O'Donnell J. E., 1994, ApJ, 422, 158 

Onken C. A., Ferrarese L., Merritt D., Peterson B. M., 
Pogge R. W., Vestergaard M., Wandel A., 2004, ApJ, 615, 
645 

Osterbrock D. E., Pogge R. W., 1985, ApJ, 297, 166 

Page D. N., Thorne K. S., 1974, ApJ, 191, 499 

Park D., Woo J.-H., Treu T., Barth A. J., Bentz M. C., 

Bennert V. N., Canalizo G., Filippenko A. V., Gates E., 

Greene J. E., Malkan M. A., Walsh J., 2012, ApJ, 747, 30 
Peterson B. M., 1993, PASP, 105, 247 
Peterson B. M., 2011, arXiv: 1109 . 4181 
Peterson B. M., et al., 2004, ApJ, 613, 682 
Pogge R. W., 2000, New Astronomy Review, 44, 381 
Pounds K. A., Warwick R. S., Culhane J. L., de Korte 

P. A. J., 1986, MNRAS, 218, 685 
Raimundo S. L, Fabian A. C., Vasudevan R. V., Gandhi 

P., Wu J., 2012, MNRAS, 419, 2529 
Richards G. T., et al., 2006, ApJS, 166, 470 
Sanders D. B., Phinney E. S., Neugebauer G., Soifer B. T., 

Matthews K., 1989, ApJ, 347, 29 
Sbarrato T., Ghisellini G., Nardini M., Taghaferri G., Fos- 

chini L., Ghirlanda G., Tavecchio F., Greiner J., Ran A., 

Gehrels N., 2012, MNRAS, 426, L91 
Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 

525 

Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 

Shen Y., et al., 2011, ApJS, 194, 45 

Shields G. A., 1978, Nature, 272, 706 

Sun W.-H., Malkan M. A., 1989, ApJ, 346, 68 

Teller R. C., Zheng W., Kriss G. A., Davidsen A. F., 2002, 

ApJ, 565, 773 
Thorne K. S., 1974, ApJ, 191, 507 
Vanden Berk D. E., et al., 2001, AJ, 122, 549 



Veron-Cetty M.-P., Joly M., Veron P., 2004, A&A, 417, 
515 

Vestergaard M., Peterson B. M., 2006, ApJ, 641, 689 

Vestergaard M., Wilkes B. J., 2001, ApJS, 134, 1 

Walter R., Fink H. H., 1993, A&A, 274, 105 

Wills B. J., Netzer H., Wills D., 1985, ApJ, 288, 94 

Woo J.-H., Urry C. M., 2002, ApJ, 581, L5 

Wright E. L., et al., 2010, AJ, 140, 1868 

York D. G., et al., 2000, AJ, 120, 1579 

Yuan W., Zhou H. Y., Komossa S., Dong X. B., Wang 

T. G., Lu H. L., Bai J. M., 2008, ApJ, 685, 801 
Zhang S. N., Cui W., Chen W., 1997, ApJ, 482, L155 
Zheng W., et al., 1995, ApJ, 444, 632 

Zhou H., Wang T., Yuan W., Lu H., Dong X., Wang J., Lu 
Y., 2006, ApJS, 166, 128 



© 2011 RAS, MNRAS OOO.[Tlf20l 



20 G. Calderone et al. 



APPENDIX B: FIGURES: SPECTRAL FITTING 

This appendix is a collection of the figures related to the 
spectral fitting procedure discussed in t )3.1l On the left pan- 
els we show the whole rest frame wavelength range, while on 
the right panels we show a detailed view on the H/3, [O ill] 
and Mgll regions. The SDSS data and associated uncer- 
tainties are shown with black squares and grey lines respec- 
tively. Also shown are the fitting models (red lines), as well 
as the individual components: the AGN continuum (black), 
the galaxy template (cyan), the jet component (as extrap- 
olated from WISE photometry, purple), the iron templates 
(orange), the broad (blue) and narrow (green) emission lines, 
and the additional emission lines (grey) . In lower part of left 
panels we show the residuals in units of data uncertainties. 
The red lines show the cumulative x?od (values on right axis) . 



APPENDIX C: FIGURES: BLACK HOLE MASS 
ESTIMATION 

This appendix is a collection of the figures related to the 
black hole mass estimation procedures described in Q 
adopting the same notation as in Fig. |4l 
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Figure Bl. Results of the spectral fitting procedure f ^S-ll Add. IB)| . 
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Figure Bl. (continued) 
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